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ABSTRACT 

We review the theoretical framework of jet photoproduction at HERA discussing 
the conceptual ideas, phenomenological models, and higher order perturbative 
calculations. Numerically, we study the physically interesting distribution of 
transverse energy within the observed jet, the real and virtual photon structure, 
and the proton structure at large x. 
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1. Theoretical Concepts 


In electron-proton scattering at HERA the dominant fraction of the scattering 
events proceeds through photons with low virtuality Q 2 .“Experimentally the electron 
is anti-tagged and remains in the beam-pipe. Theoretically the lepton tensor and 
phase space can then be factorized with the Weizsacker-Williams approximation!, 
which was improved recently through power-suppressed terms by Frixione et al.! The 
virtuality of the photon is less than 4 GeV 2 , and it retains a fraction y E [0.2; 0.85] 
of the incident electron energy E e = 27.5 GeV. The proton energy is 820 GeV. 

Jet production in photon-proton collisions was first calculated in leading order 
(LO) of perturbative QCD by Owens!. The two-fold nature of direct and resolved 
processes and their separation on a kinematical basis were pointed out as a means to 
test the underlying partonic dynamics. It is thus possible to study the distribution 
of partons in the initial photon and proton and the spin of the exchanged particle. 

Since then much effort has been spent on improving the tree-level understanding. 
Phenomenological models have been implemented into the Monte Carlo event genera¬ 
tors PYTHIA0, HERWld, and PHOJETl employing parton showers, fragmentation 
models, and multiple interactions. Next-to-leading order (NLO) calculations for inclu¬ 
sive single-jet! and dijet production! provide perturbative correction factors, reduce 
the scheme and scale dependences, and allow for an implementation of jet definitions. 
Recently the first NLO analysis of the transition from real to virtual photoproduction 
of jets has been published!. 

Although several parametrizations exist for the parton distributions in the photon 
in NLO0, the gluon is still poorly constrained by the data from 77 scattering at e + e~ 
colliders and needs complcmentarvmeasurements from 7 p. The situation is better for 
the parton densities in the protonlHl, where photoproduction data might improve the 
data from deep inelastic and pp scattering at intermediate to large values of x. Except 


“From 1994-1997, HERA operated with positrons instead of electrons. Since we will only be con¬ 
cerned with neutral current exchange, we use the term “electron” for positrons as well. 



where indicated we use the GS96 and CTEQ4M parametrizations for the photon and 
proton, respectively. 

2. Phenomenological Models 

Perturbative QCD describes only the hard partonic scattering process and the 
scale evolution of the hadronic structure functions. The link between colored partons 
and real hadrons belongs to the non-perturbative domain and has to be fitted to data 
or guessed from phenomenological models. These models are implemented in Monte 
Carlo generators and compared in Table 1. 


Table 1. Properties of different Monte Carlo generators. 


Monte Carlo 
Generator 

Parton 

Showers 

Fragmen¬ 

tation 

Multiple 

Interactions 

PYTHIA 

Initial+Final 

String 

Hard (optional) 

HERWIG 

Initial+Final 

Cluster 

Soft (optional) 

PHOJET 

Final 

String 

Soft and Hard 


As a first step, one can attach additional angularly ordered partons to the hard 
process until the original parton reaches a maximum virtuality of Q„ iax < E f>. In 
addition and deviating from the parton model, intrinsic transverse momenta are al¬ 
lowed for the partons in the hadrons up to kr < /c™ ax thereby introducing a second 
phenomenological parameter. 

The second step consists in the fragmentation of the proliferated partons into 
hadrons. The Lund string model confines the color field between quarks into a color 
flux tube. The energy increases proportional to the distance between the quarks until 
the string breaks up and new quark-antiquark pairs are created. Hadrons are formed 
when the energy is too low for the string to break up according to a fragmentation 
function with two free parameters. Gluons appear as excitations and produce kinks in 
the string. Alternatively the cluster model starts with the splitting of gluons into qq 
pairs and the subsequent formation of color neutral clusters. Heavy clusters cascade 
into light clusters which then transform isotropically into hadrons. 

Finally multiple interactions seem to be important at HERA to describe transverse 
energy flow and cross sections in the direction of the proton remnant at low Et■ They 
increase the multiplicity and energy flux and model interactions between the photon 
and the proton remnant. Secondary interactions are softer than the first scattering 
defined by Et and can be of partonic or soft nature. 

3. Next-to-Leading Order Calculations 


The undesirable drawback of phenomenological models is the large number of free 
parameters that have to be tuned to data. These are not present in next-to-leading 




order calculations where one calculates virtual corrections with internal particle loops 
and real corrections with soft and collinear radiation one order higher in the strong 
coupling constant a s . The singularities can then be controlled through dimensional 
regularization and removed consistently through renormalization and factorization 
procedures. This avoids unphysical cut-off parameters, reduces the scheme and scale 
dependence, and allows for an implementation of jet definitions. 

Real corrections are most easily calculated with the phase space slicing method. 
After approximation of the invariants and partial fractioning of the 2 —> 3 matrix 
elements, we factorized the Born process and integrated the remaining singular ker¬ 
nels analytically up to an invariant mass cut-off y cut . Aurenche et al.! employed a 
transverse energy cut-off for the soft/collinear initial state, which does not appear in 
a related method for the single-jet case, and a cone cut-off for the collinear final state? 
Harris and Owens integrated the soft and collinear regions separately!. The real emis¬ 
sion outside the cut-offs is integrated numerically. This removes the dependence on 
the technical cut-off and introduces the experimental jet definition. 

The subtraction method relies on a point-by-point subtraction of singularities in 
the numerical integration and has only been applied to the direct process by BodekeiS. 
The resolved process could, however, be adapted from an existing program in pp 
scattering^. 

4. Jet Definition Uncertainties 


In hadronic collisions cluster algorithms of the JADE type combine particles not 
only into the hard jets but also into the remnant jets. This is avoided if a cone of 
size R = 1 in azimuth-rapidity space is used to define a jet. The Snowmasdi! accord 
determines the jet axis from the Et weighted directions of all particles in the cone. 
Unfortunately it contains a number of ambiguities which render a matching between 
theory and experiment difficultly. 

Contrary to fixed cone algorithms, iterative cone algorithms can merge overlapping 
jets. This is not described by a NLO calculation with just three final state partons. 
If two partons have a distance between R and 2R, they can be counted as one or two 
jets, and one has to avoid double-counting. The phenomenological parameter R sep 
defining the distance of two partons can be used to model the narrower jets found by 
iterative cone algorithms. 

In a recent studyEl, we used this R sep parameter to describe jet shapes p(r ) as 
measured by the ZEUS, CDF, and DO collaborations. The result is compared to pre¬ 
liminary ZEUS data from 19940 which was obtained with an iterative cone algorithm 
and is shown in Figs. 1 and 2. It depends moderately on the rapidity ij of the jet 
(Fig. 1) being broader in the forward direction and narrows slightly with the trans¬ 
verse energy Et of the jet (Fig. 2). Whereas the curves with R sep = 2 correspond to 

b See also the contribution by M. Fontannaz in these proceedings. 



no R sep and predict too broad jets, an average value of -R sep = 1.4 describes the data 
rather well. 

The uncertainties of the Snowmass cone definition are circumvented in the hadronic 
implementation of the Er-cluster algorithm^ with the same EV-weighted recombina¬ 
tion scheme as for Snowmass, but a different jet condition. One considers the distance 
of two particles in rj —■ <f> space. This corresponds to a unique value of R sep = R in 
theory and experimentally assigns every hadron to a unique jet. 

5. Real and Virtual Photon Structure 

The determination of the structure of the photon is clearly one of the most im¬ 
portant physics goals in photoproduction. Since the cross section drops rapidly with 
the transverse energy, most events at HERA have been observed at small Et so far. 
They are dominantly produced by resolved photons and give access to the poorly 
constrained small-x region and the gluon density in the photon. 

However at low Et a separation of hard and soft physics is experimentally and the¬ 
oretically difficult. Hadronization effects, jet definition uncertainties, energy pedestals 
around the jets, and multiple interactions between the remnant jets may play an 
important role. This can be seen in Figs. 3 and 4, where we compare our NLO 
calculation to preliminary 1994 data from ZEUS, again obtained with the iterative 
cone algorithmEJ. At Et > 14 GeV, there is a clear excess of data over theory in 
the forward direction, which decreases continuously for larger Et- We do not show 
a comparison with data obtained with a smaller cone size R — 0.7 which shows no 
excess in the forward region. Obviously the hard jets are better separated from the 
underlying event for smaller cones and there is no need for multiple interaction effects. 

A second conclusion from Figs. 3 and 4 is that the jet definition uncertainty is of 
comparable size to the photon structure function uncertainty. Thus one either has to 
rely on the fitted value of R sep = 1.4 from jet shapes or use the hr cluster algorithm. 

In the backward region, where direct and quark initiated processes dominate, there 
is fairly good agreement with the data with a slight tendency of GRV to overestimate 
the measurement. This can be understood from the bigger quark contribution of GRV 
compared to GS96 at large values of x. 

The kinematics of the partonic subprocess can be better constrained in dijet cross 
sections than in single cross sections, e.g. with the help of the variable x® BS = 
+ Er 2 e~ m ) / {2yE e ). As long as one integrates over different bins, this vari¬ 
able is infrared-safe and there is no need to abandon it as put forward by Aurenche 
et al.l With the help of this variable, ZEUS separated direct and resolved regions 
experimentally and studied the dependence of the cross section on the average ra¬ 
pidity of the observed jets@. Unfortunately they constrained the Et of both jets to 
the same minimal values. The theoretical predictions are then not infrared safe and 
depend on a phenomenological (Klasen and Kramer) or the technical cut-off (Harris 


and Owens)!. This analysis employed the kr cluster algorithm with R = 1 and ex¬ 
hibited again the excess of data over theory for the resolved dominated regions at low 


and E t . 


A new preliminary analysis of 1995 ZEUS data has been presented recentl 
They measured the symmetrized dijet cross section da / dE T drjidr}2 with the leading 
Et > 14 GeV and E T2 >11 GeV in different rapidity bins and the full and upper 
region of x° BS . The result is infrared safe and compared to our NLO predictions in 
Fig. 5 as a function of Et and in Fig. 6 as a function of 772 . The general agreement is 
good for both GRV and GS96, even when both jets are in the forward region and the 
complete range in x® BS is covered. This may be due to the fairly large cuts on the 
Et of both jets, which suppresses the underlying event. In the backward region, the 
calculation lies above the data and GS96 is slightly favored as was the case for the 
single-jet cross sections. In addition the systematic errors not shown here are quite 
prominent in the backward region. Since the data using the kr algorithm were not 
yet available, we simulated the iterative cone algorithm with the optimized value of 
R sep — 1.4 for R — 1. This value should be lower in the backward region leading to 
a reduced theoretical prediction and better agreement with the data. 

A new NLO program for virtual photoproductionl! allows us to study the the tran¬ 
sition to photons with virtuality P 2 . Apart from using the unintegrated Weizsacker- 
Williams approximation, the main difference consists in the analytic integration of 
the virtual photon initial state singularity. This singularity is then factorized and 
leads to a scheme- and P 2 -dependent finite contribution which reduces to the real 
expression as P 2 —» 0. How this leading-logarithmic singularity compares to the two 
existing up-quark parametrizations in the virtual pliotoriEU is shown in Fig. 7. At low 
P 2 = 1 GeV 2 the evolution up to 50 GeV 2 produces many more up-quarks than just 
the single pair predicted perturbatively. At large P 2 = 10 GeV 2 , GRS and SaSIM 
agree with the perturbative box-diagrami!. 

ZEUS have also published P 2 dependent data on the ratio of resolved over direct 
contributions disentangled with the x° BS cut at 0.7^1. As can be seen in Fig. 8 the 
NLO effects are large at Et 1 ,Et 2 > 4 GeV. The lowest P 2 point for real photopro¬ 
duction has much better statistics than the virtual photoproduction data and lies 
above the theory. This is to be expected since we are again in the soft region where 
remnant-remnant interactions are assumed to be important. At larger P 2 , theory and 
data agree well in shape and normalization. We use the SaSIM LO parametrization, 
since NLO parton densities are not available in parametrized form and GRS does not 
contain charm. 


6. Proton Structure at large x 

The excess of liigh-Q 2 events in DIS at HERA have triggered a lot of speculation 
about new physics such as leptoquarks, P-parity violating supersymmetry, or contact 


interactions. This kinematic regime can and should also be tested in photoproduction 
at large E T , e.g. by measurements of da/dE T , da/dE^, or da/dMjj. An unambiguous 
determination of new physics requires a precise knowledge of the Standard Model 
background. Due to phase space limitations, NLO corrections become increasingly 
important at the boundary of phase space at large E Tl whereas the uncertainty from 
the photon structure function becomes negligible. This can be extracted from Fig. 9, 
where the fractional contribution of direct photon-quark scattering accounts for more 
than 90% of the total single-jet cross section above E? = 100 GeV. Therefore a precise 
knowledge of the quark distribution in the proton is required for x > 0.4. 

In Fig. 10, we demonstrate the increase in the cross section that can be obtained 
with the larger quark densities of MRS(J’) which were designed to describe the excess 
of CDF pp data at large eM- A factor of 1.5 seems possible here. Larger gluon 
distributions as proposed by CTEQ4(HJ)0 only give a rise of 15%. It should be 
mentioned that MRS(J’) fail to fit the low energy BCDMS data, whereas CTEQ4(HJ) 
are still consistent. Finally, a study of the scheme dependence might indicate further 
uncertainties here or in the gluon density. 

7. Conclusions 

A wealth of new and precise data coming from HI and ZEUS in photoproduction 
has triggered increased theoretical interest. Three independent NLO calculations 
for real dijet production and the first NLO program for virtual photoproduction are 
now available. The interpretation of jets is improved with smaller cone sizes and 
the kx cluster algorithm. This will allow for first stringent tests of the real photon 
parton density, especially of the gluon. More data and a parametrization of NLO 
parton densities are needed in virtual photoproduction. The access to large transverse 
momenta offers interesting studies of the proton structure at large x. 
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Fig. 1 . Jet shape p(r) for single-jet photoproduction integrated over Et >14 GeV and four different 
regions of r : ]. We compare our results using the Snowmass convention with R= 1 and three different 
values of Rse P to preliminary 199j data from ZEUS. 
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Fig. 2. Jet shape p{r) for single-jet photoproduction integrated over —1 < 77 < 2 and four different 
regions of Et- We compare our results using the Snowmass convention with R = 1 and different 
values of R sep to preliminary 199j data from ZEUS. 




















Fig. 3. Rapidity dependence of the single-jet photoproduction cross section integrated above Et > 14 
GeV. We compare our NLO prediction with GRV and GS96 photon parton densities and the two 
extreme values of R sep = 1R,2R to preliminary 1994 data from ZEUS. 



Fig. 4. Same as Fig. 3 for Et > 21 GeV. 
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Fig. 5. Et dependence of the symmetrized dijet photoproduction cross section integrated over different 
rapidity bins. We compare our NLO prediction with GRV and GS96 photon parton densities and 
the full and upper range of x® BS to preliminary 1995 data from ZEUS. 
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Fig. 7. Up-quark distributions in the virtual photon as function of x. We compare the GRS and 
SaSIM parametrizations with the perturbative result for P 2 = 1 and 10 GeV 2 l£3. 



Fig. 8. Ratio of resolved over direct contributions defined with the x° BS cut at 0.15 for E , Et 2 > 4 
GeV using the SaSIM virtual photon structure with four flavors. 


























FractiorK.it Difference tg T 3 Fractional Contribution 















